Optically controlled spin-glasses in multi-qubit cavity systems 
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Recent advances in nanostructure fabrication and optical control, suggest that it will soon be 
possible to prepare collections of interacting two-level systems (i.e. qubits) within an optical cav- 
ity. Here we show theoretically that such systems could exhibit novel phase transition phenomena 
involving spin-glass phases. By contrast with traditional realizations using magnetic solids, these 
phase transition phenomena are associated with both matter and radiation subsystems. Moreover 
the various phase transitions should be tunable simply by varying the matter-radiation coupling 
strength. 

PACS numbers: 42.50.Fx, 75.10.Nr, 32.80.-t 



Condensed matter physicists are keen to identify ex- 
perimental realizations of Ising-like Hamiltonians in- 
volving populations of interacting two-level objects (i.e. 
spins) |l| . Exotic phases such as a spin-glass are of par- 
ticular interest [21 . Experimental studies of spin-glasses 
have focused on solids containing arrays of magnetic ions 
[ll. However the laws of Nature limit the range of ex- 
otic behaviors that such solids can exhibit, since it is 
very hard to engineer the magnitude, anisotropy, range 
and/or disorder of the spin-spin interaction in such sys- 
tems. In the seemingly unrelated fields of atomic, nanos- 
tructure and optical physics, there have been rapid ad- 
vances in the fabrication and manipulation of effective 
two-level systems (more commonly referred to as qubits) 
using atoms, semiconductor quantum dots, and super- 
conducting nanostructures y, la la Li la Is! ■ In particular, 
controlled qubit-cavity coupling has been demonstrated 
experimentally between such systems and a surrounding 
optical cavity 0, |a la S la ISl lldl- ^^ quantum dots 
systems, in particular, the effective spin-spin (i.e. qubit- 
qubit) interaction can in principle be tailored by adjust- 
ing the quantum dots' size, shape, separation, orienta- 
tion and the background electrostatic screening. Further- 
more, the interaction's anisotropy can be engineered by 
choosing asymmetric dot shapes. Disorder in the qubit- 
qubit interactions will arise naturally for self-assembled 
dots, or can be introduced artificially by varying the in- 
dividual dot positions during fabrication [sl, |llj . 

Motivated by these recent experimental advances, we 
study the phase transitions which could arise in such 
multi-qubit-cavity systems. We uncover novel realiza- 
tions of spin-glasses |3| in which the phase transition phe- 
nomena are associated with the spin (i.e. matter) and 
boson (i.e. photon) subsystems. The resulting phase di- 
agrams can be explored experimentally by varying the 
qubit-cavity coupling strength A, e.g. by re-positioning 
the center of the cavity around the nanostructure array, 
changing its orientation, or tuning the cavity Q-factor 
|2J. In addition to opening up the study of these im- 
portant condensed matter systems to the nano-optical 
community, our results help strengthen the theoretical 



connection which seems to be emerging between multi- 
qubit-cavity systems and spin-spin systems |l2l |l3J . 




FIG. 1: (color online) Schematic diagram showing a possi- 
ble optical realization of a spin-glass system within a multi- 
qubit-cavity system comprising elongated quantum dots. The 
quantum dots are orientated along the a;-axis and have non- 
uniform spacings, giving rise to disorder in the qubit-qubit 
interaction {Jij}. 

Given the above experimental considerations, we will 
int rodijce a gfineralization of the well-known Dicke model 
|l2. Ha. I14L Tl5l | in order to describe N interacting two- 
level systems in a cavity field: 
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where the operators a, a^ and a ■ , a^ correspond to the 
photon field and quantum dot j respectively. The inter- 



action term Hi 



T,i<j Jij'^i'^j or Hi 



Z-ii<i "^ij 



(or equivalently J2t<j Jij^l^])- 
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An interaction in the 
direction is present naturally 4,'. In the case of the 



x-directed interaction, each quantum dot can be engi- 
neered to have an elongated form along the x-direction, 
by biasing the growth process along this direction. Ap- 
plying an electric field along x, will then create large 
dipole moments in that direction. One can use undoped 
dots, in which case the dipole results from the exci- 
ton, or doped dots, in which case the dipole originates 
from the conduction-subband electron biased along x. A 
schematic plot of a possible realization is depicted in Fig. 
1. Kiraz et al. [a|, for example, have already built physi- 
cal realizations of such multi-qubit-cavity nanostructure 
systems. In particular, Ref. [Ej shows a scanning electron 
micrograph of the GaAs microdisk nanostructure system 
comprising a disordered array of quantum dots embed- 
ded within an optical cavity. The cavity has a diameter of 
4.5^TO and produced Q-factor values exceeding 1800 [5j- 
The cavity was built with a collection of InAs quantum 
dots at irregular locations fixed during the growth, and 
hence necessarily features disordered dot-dot couplings 
{Jy}. Another suitable experimental set-up has been 
provided by Imamoglu et al. [a| in which a quantized 
cavity mode and applied laser fields are used to mediate 
the interaction between spins of distant, doped quantum 
dots. This leads to cavity-assisted spin-flip Raman tran- 
sitions and hence pairs of quantum dots can be coupled 
via virtual photons in the common vacuum cavity mode 



The dot-dot coupling terms { Jij } will have an inherent 
disorder in self-assembled dots - alternatively, such disor- 
der can be built in during growth by varying the dot-dot 
separations. We make the reasonable assumption ^^ that 
the disorder will be Gaussian: 



P{J.j) 



( ^ 



V27rJ2 



exp 



4^ 

' 2J2 



First we consider Hi ~ J2i<i Jij^t'^^ji noting that the 
same method can then be used to solve for spin-spin in- 
teractions in the other directions. To obtain the ther- 
modynamical properties of the system, we introduce 
the Glauber coherent states \a) of the field [l6j where 
a|a) = a\a), {a\a^ = {a\a* . The coherent states are 
complete, / ^^\a){a\ = 1. In this basis, we may write 
the canonical partition function as: 

Z(Ar,T) = V/^(s|(a|e-^«l 



as. 



As in Ref. [l5J, we adopt the following assumptions: 

1. a/VN and a' /vN exist as TV ^ oo; 

2. limTv^oo linifl^oo Y.r=o ,?"'' can be inter- 
changed. 

We then obtain: 
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Performing the Gaussian integration of Im(Q:) and setting 
w = Re (a) yields the partition function: 
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We set Jij — > —Jij, e = — 2/i and ^ = ^^ for convenience 
and rescale the photon field a^ -^ ^/Nb^ such that w -^ 
VnW. We then use the Trotter- Suzuki method Tr|, and 
obtain the free energy from the replica trick: 



-/3F = [ln{Z)] = lim 



[Z"] - 1 



n^O n 

where [Z"] is the configurational average of the n repli- 
cated partition functions. Rewriting the trace and using 
the method of steepest descents, we can proceed in the 
usual manner for spin-glass problems. In the thermody- 
namic limit N —f oo, with the limit n — > with N kept 
large but finite, wc obtain the free energy, f — jj, as: 
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with a being the replica index and t = 1, 2, 3, • • • , P the 
label for the Trotter direction. We then substitute for 
the following stationary values: 
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The replica symmetric assumption [l7| yields x'tv — Xtt' , 
q"^, = qtt', W" = W, and with the static approximation 
we then obtain xw = X and qtt' = q- The validity of this 
approximation has been discussed at length in Ref. |l7| . 
Hence: 
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y^Jexp J — ^ L We then calculate the 

and perform power series expansions in the limit 
P -^ oo. Finally, we expand out the last term in the 
free energy, take the limit n — > and replace /z and F by 
their original values fi = —^ and F — — 2A to obtain: 
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with 

h ^ pj [^zi + Vx - qz2 - - j , 

having absorbed the factors ^ and 4 into e and A respec- 
tively. The requirement that the free energy is stationary 
yields the following self-consistent equations: 
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Adding AiJ = haf to the Hamiltonian such that e 
(e -I- h) and taking the limit /i — > 0, yields: 



Pdh 
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An intricate array of phase transitions arises in both 
the matter and the optical sub-systems. Furthermore, 
these transitions in the matter and optical sub-systems 
are inter-dependent. Here we summarize the main re- 
sults for a general interaction Hi — X]i<7- Jij'^f'^j + 
J2i<j '^ij'^i'^j- If tl^s interaction is only in the z- 
direction, a ferromagnetic-paramagnetic phase transition 
emerges in the matter sub-system - but there is no spin- 
glass phase anywhere in the phase diagram. This is to be 




FIG. 2: A plot of the magnitude of the super-radiant order 
parameter, W , against A with /3 = 2, e = 1 and J = 2, and 
a spin-spin interaction in the z-direction. The gap between 
the two non-zero W phases disappears as we decrease the 
temperature, T. 



expected, as at any finite value of e there is a non-zero 
contribution to the magnetisation as e acts as a longitu- 
dinal term to the Ising interaction 2J. As such, we would 
only expect a spin glass phase in the limit of e ^ 0. How- 
ever, novel behavior emerges in the optical sub-system in 
the form of two sub-superradiant phase transitions, as 
shown in Fig. 2. This effect disappears as we increase 
/?. We see that the system goes through a first order and 
then a second order phase transition as we increase A |l,ll | . 
In particular, the sharp nature of both transitions and 
the narrowing of the gap between them as (3 increases, 
suggests that an accurate low-temperature thermometer 
could be constructed based on this effect. 

The addition of the x-directional interaction to the 
original Dicke Hamiltonian yields a number of interest- 
ing results. In particular, the condition for the system 
to exhibit a sub-superradiant phase transition is drasti- 
cally changed from the case without interactions. If J 
is large, the matter sub-system swamps the optical sub- 
system and the system is unable to reach superradiance. 
Consequently the magnetization is always zero. However, 
lowering J yields a superradiant phase - and by varying 
e wc uncover a spin-glass phase that moves into a para- 
magnetic phase at higher e. Most remarkably, there are 
two different phase diagrams for the system depending on 
the value of the qubit-cavity coupling A at a given J. The 
corresponding phase boundaries are indicated in Fig. 3. 
The complex behavior in the region between these two 




FIG. 3: (color online) Schematic plot of the phase diagram 
of our system when a spin-spin coupling is present in the x- 
direction. As we increase A, the phase diagram changes from 
one which resembles that of an Ising spin-glass model with a 
transverse field, to one with only ferromagnetic and param- 
agnetic phases (i.e. the spin-glass phase has disappeared). 



limits will be discussed in depth elsewhere. 

To summarize, our work shows that novel phase tran- 
sition phenomena will arise in suitable optical realiza- 
tions of generalized spin-spin Hamiltonians. In addition 
to their intrinsic theoretical interest, we hope that our 
findings might motivate experimental efforts toward ex- 
ploring these predicted phases. 
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